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ABSTRACT 

We used the Spitzer Infrared Spectrograph (IRS) to study the mid-IR properties and investigate 
the energetics of 24 Submillimctcr Galaxies (SMGs). This is the largest sample of SMGs observed 
with Spitzer IRS and covers the full extent of properties of the radio-identified population of SMGs 
in the redshift range of z ~ 0.65-3.2. We observe broad emission features from Polycyclic Aromatic 
Hydrocarbons (PAHs) in more than 80% of our sample. We find that the median mid-IR spectrum is 
well described by a starburst component with an additional power law, F v ~ v~ 2 , likely representing a 
< 32% AGN contribution to the bolometric luminosity. Our results thus confirm that starburst activity 
dominates the bolometric luminosity in SMGs. We find that SMGs show weaker silicate absorption 
at ~ 9.7/im than local ULIRGs. We also find stronger 6.2-//m PAH emission in SMGs (relative to 
the 7.7/im PAH feature) than in local nuclear starbursts, which may be attributed to lower extinction 
by ice along the line of sight to SMGs. This suggests that the continuum and PAH emitting regions 
of SMGs are less obscured than in local starbursts and similarly luminous low-redshift ULIRGs. We 
interpret these results as evidence for a more extended distribution of cool and warm dust in SMGs 
compared to the more compact emitting regions in local ULIRGs and starbursts. Together these 
results suggest that SMGs are not simple high-redshift analogs of local ULIRGs or nuclear starbursts, 
but instead they appear to have star formation which resembles that seen in less-extreme star-forming 
environments at z ~ - suggesting their intense activity is distributed across a far larger region than 
the ~ 1-kpc nuclear bursts in local ULIRGs. 

Subject headings: infrared: galaxies — galaxies: starburst — galaxies: AGN — galaxies: submillimeter 



1. INTRODUCTION 

It is over a decade since the discovery of a pop- 
ulation of high-redshift galaxies identifi ed through 
their submillimeter (submm) emission (ISmail et al.l 
19971 iBareer Cowie feijanderj p9l lEales et al.l ]1999j 
Bertoldi et al.ll2000t ICowie et al . 2002: S cott et alJ | 2002| 
Borvs et al.l 120031: 1 Webb et al.ll2003t ICoppin et alj|2005l 
Younger et al.l |2007( ). Their submm selection sug- 
gests that these are strongly star-bursting systems. 
With implied star-formation rates (SFRs) of ~ 100- 
1000M Q yr _1 , these dust-enshrouded submm galaxies 
(hereafter, SMGs) may make a sign ificant contribution to 
the global SFR density at z ~ 2-3 (|Chapman et al.ll2005l 
hereafter C05). The absorbing dust that makes them 
such prodigious submm emitters also makes them quite 
optically faint and renders follow-up studies at shorter 
wavelengths challenging. The study of SMGs has been 
facilitated by the detection of a large fraction of them 
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as //Jy radio sources (jlvison et al.l [2002) and more re- 
cently through millimeter-wave intcrfcromctry. This has 
allowed the subsequent measurement of their spectro- 
scopic redshifts for fairly large samples (C05). 

With a mean redshift of < z >^ 2.2, the redshift dis- 
tribution of the radio-identified SMGs in C05 coincides 
with the global peak epoch of quasar activity. The con- 
nection between SMGs and active galactic nuclei (AGNs) 
has been probed in recen t near-IR and X-ray s tudies. 
Near-IR spectroscopy by ISwinbank et al.l (|2004h show 
that broad Ha lines (FWHM > 1000 kms" 1 ) are often 
present in these galaxies, while deep X-ray studies using 
the sensitive Chandra Deep Field-North (CDF-N) Sur- 
vey suggest t hat at least ~ 28-50% of SMGs host an ob- 
scured AGN (| Alexander et al.ll2005h . A number of SMGs 
that display no AGN s ignatures in the rest-frame optical 
(jSwinbank et al.ll200l were classified as AGNs based on 
X-ray observations ( Alexander et al.ll2005h . This is likely 
due to geometrical effects in which the broad-line region 
of the AGN remains hidden in the optical by intervening 
obscuring material. At high X-ray energies, direct X-ray 
emission is detectable through even very high column 
densities, allowing the direc t detection of AGN at very 
high obscurations. However. I Alexander et al.1 (|2005l ) find 
that the majority of AGNs in SMGs are more modestly 
obscured, with column densities of N# > 10 23 cm~ 2 . A 
similar mix of AGN and starburst activity is seen in lo- 
cal ultra-luminous infrared galaxies ( ULIRGs, with to- 
tal IR luminosities L m > 10 12 L Q ; ISoifer et al.l 119871 : 
iSanders et al.l fl988). most of which have b een shown to 
be co mposite AGN-starburst systems fe.g.. lArmus et all 
I2007T ). The SMGs have IR luminosities which are com- 
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parable to ULIRGs at low redshift, prompting the ques- 
tion as to whether SMGs are high-redshift analogs of 
ULIRGs and hence whether we can learn about the phys- 
ical processes within SMGs from studies of local ULIRGs 
(|Tacconi et all 12001 ) . 

The potential presence of luminous Compton-thick 
AGNs (N H > 10 24 cm" 2 ) in those SMGs with no X-ray 
AGN signature remains a significant caveat to these re- 
sults. Furthermore, the samples of SMGs with the neces- 
sary ultra-deep X-ray observations to reveal the presence 
of highly obscured AGN are still small. Hence it is possi- 
ble that a small fraction of luminous, but C ompton thick, 
AGN lurk within the SMG population (see lCoppin et all 
l200l) . 

Rest-frame optical emission provides direct insight into 
the stellar emission and ionized gas of a galaxy, but 
suffers from obscuration due to intervening dust. At 
longer wavelengths, the mid-IR emission arising from 
the dust itself provides an indirect insight into the dust- 
enshrouded nature of SMGs and suffers from much less 
obscuration than the shorter, optical wavelengths. We 
therefore designed a program to follow up a large sam- 
ple of 24 SMGs in the mid-IR using Spitzer's Infrared 
Spectrograph (IRS: lHouck et al.ll2004l ) in an effort to an- 
swer the following questions for a large and representa- 
tive sample of the SMG population: Are SMGs compos- 
ite AGN-starburst systems? To what extent does AGN 
activity contribute to the total infrared output of these 
galaxies? Is there a spectrum of varying levels of AGN 
activity across the population? Are local ULIRGs good 
analogs for the mid-IR emission of SMGs? 

The main components contributing to the mid-IR spec- 
trum of a galaxy are: thermal dust continuum, emission 
from vibrational modes in Polycyclic Aromatic Hydro- 
carbons (hereafter, PAHs) and other atomic and molecu- 
lar lines. The continuum emission at longer mid-IR wave- 



lengths, A > 12/im, arises from emission by very small 
grains of dust (VSG; < lOnm) found around obscured 
AGN or star-forming regions. This is often referred to 
as the VSG continuum, or the warm dust continuum 
O^Dust < 250 K). At shorter wavelengths, A < 6^m, the 
continuum traces emission from dust heated to signifi- 
cantly hotter temperatures (Trj ust > 500 K) likely due 
to its close location near an AGN or possibly a hot, nu- 
clear starburst region. This is what we refer to as the 
hot dust continuum. PAH molecules (< few nm) consist 
of chained benzene rings, associated hydrogen and other 
trace elements, such as Si and Mg. The line emission in 
the mid-IR waveband is dominated by PAH molecules 
which are excited by the UV photons that are copi- 
ously available in star-forming regions. The main PAH 
emission features arise from the bending and stretching 
of skeletal C-C or peripheral C-H bonds and are ob- 
serve d at rest- frame 6.2 , 7.7, 8.6, 11.3, 12.7 and 17 fim 
(e.g.. lDraine fc Lil l20071 It has been shown locally that 
stronger PAH feature s are associa ted with regions of in- 
tense star formation fHelou 199j|). 

With the unprecedented sensitivity of Spitzer IRS, it 
has been possible to explore the mid-IR region of galaxies 
at high re dshift down to continuum levels of S 24 am > 0-1 
mJv (e.g.. lYan et al.l[2005l. [20071: ILutz et al.ll2005l 120071: 
IDesai et all 120061: jRigopoulou et all 120061 : iTeplitz et all 
l2007t ISiana et all 12008). At lower redshifts, investiga- 
tion of the mid-IR properties of local galaxies with IRS 
provides for a sample of detailed templates against which 
high- r edshift sources can be com p ared ( e.g.. iSpoon et al. 
200 l lATmus et al l 12001 [200l [20071: IWeedman et al. 
2005 : iBrandl et alj|200a ISturm et alj[2006t IDesai et al. 
20071 ). A wide range of mid-IR spectra has been uncov- 
ered for high-redshift sources, ranging from continuum- 
dominated spectra wit h no PAH features to PAH- 
dominated spectra (e.g.. IWeedman fc Houckl l2008f). Be- 
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Fig. 1. — (Left:) The distribution of optical redshifts and IR luminosities (8-1000/im) of the parent sample of radio-identified SMGs 
compiled by C05 with our IRS sub sample of SMGs h ighlighted. We use IR luminosities based on the values reported by C05, corrected 
down by a factor of two, following Kovacs ct al. (2006). We plot a representative error on the IR luminosities, corresponding to ~ 20% 
(C05). (Right:) The distribution of radio and submillimeter fluxes for the full C05 sample and those in our subsample observed with 
IRS. The histograms show the distribution in parameters of the full sample (dotted line) and our IRS subsample (solid line). These plots 
demonstrate that the IRS sample is representative of the parent sample of radio-identified SMGs, covering the full range in radio and 
submm properties, redshift and IR-luminosity. 



tween these two extremes there is a myriad of composite 
spectra displaying a sign ificant continuum with super- 
posed PAH features fe.g. JYan et al.ll2007l) . 

We presented the Spitzer IRS results for the 
first five SMGs observed i n the program in 
iMenendez-Delmestre et all (|2007T ). which complemented 
ea rlier IRS re s ults f or two SMGs at z ~ 2.8 presented 
bv lLutz et all (I2005D. Since then, IValiante et all (|2007l ) 



and iPone et alj ( 2008f ), obtained IRS spectra of nine 
SMGs in blank field and cluster lens surveys and of 13 
SMGs from the GOODS-North Field, respectively. In 
this paper we present the results for our full observing 
program of 24 SMGs, comprising the largest sample 
of SMGs observed to date with Spitzer and nearly 
doubling the number of SMGs observed with IRS. We 
discuss our sample selection and observing strategy in 
§2] The steps comprising the reduction and analysis of 
the spectra are discussed in §21 and SJ21 Our results are 
presented in §1] and discussed in detail in £[5] We give 
our conclusions in £j6j We assume a ACDM cosmology, 
with H = 71kms- 1 Mpc" 1 , Q M = 0.27 and n A = 0.73. 

2. OBSERVATIONS AND REDUCTION 

The SMGs in our program were selected from the sam- 
ple of 73 spectroscopically-confirmed, radio-identified 
SMGs in C05. We aimed for full coverage of the range in 
redshifts, radio, submm and total IR luminosities present 
in the C05 sample (see Fig.[T]) to build a sample represen- 
tative of the radio-identified SMG population in general. 
The SMGs in our sample have optical redshifts 0.6 < z < 
3.6, radio and submm fluxes 40/iJy< Si.4GHz ^ 240/iJy 
and 4 mjy< Ssso^m ^ 12 mJy (C05), and IR luminosities 
2 x 10 11 L Q < -Ls-iooo^m < 6 x 10 13 L Q , where the latter 
are based on the bolometric values reporte d by C05, cor- 
rected down by a factor of two, following iKovacs et al.l 
(2006). We note that several of the SMGs from the C05 
spectroscopic sample fall in the GOODS-North region 



and were included in the IRS studies focused on that 
field (jPope et al.ll2008l ). Together with our sample, these 
are effectively part of the larger IRS study of the C05 
survey which covers a total of 33 SMGs (or 45% of the 
parent C05 sample). 

While we aimed to span the full range in multiwave- 
length properties of radio-identifed SMGs, to have a re- 
alistic chance of detecting the galaxies with IRS we se- 
lected the 24 SMGs with higher estimated 24/im fluxes. 9 
As Spitzer 24/im observations were not available at the 
time of the proposal, we estimated 24/tm fluxes for the 
SMGs in C05 from fitting the spectral energy distribu- 
tions (SEDs) to the radio (1.4 GHz), submm (850/xm) 
and available optical photometric data points (see C05 
and references therein). We assumed a power-law fit, 
S v ~ v~ a to the mid-IR with a = —1.7. We note that 
selection of 24/im-bright sources (albeit estimated) may 
introduce a two-sided selection effect. Selection based 
on bright 24/im-fluxes will preferentially pick out the hot 
dust emission in AGNs. However, at redshifts z ~ 2, the 
7.7/im-PAH feature is rcdshifted onto the 24/im-band: 
this is a prominent emission feature, associated with star 
formation activity and will result in a boosted 24/im 
flux. We discuss this potential selection bias in more 
detail in §4.5.21 From a parallel mid-IR imaging sur- 
vey of the radio-identified SMGs i n C05 with the Multi- 
band Imager for Spitzer (MIPS) bv lHainline et all (|2009l ) 
our IRS subsample has MIPS 24 /im band fluxes in the 
range 5*24 Mm ~ 0.09-0.85 mJy with a median value of 
< S 2 4 M m >= 0.33 ± 0.18 mJy. 

The details of our observations are summarized in 
Table [TJ We observed all our targets using the low- 
resolution (i? ~ 57-126) Long-Low (LL) observing mod- 
ule of IRS. The slit width for this module is ~ 10.7", 

9 We used the online tool SPEC-PET to calc ulate S/N estimates 
( http: / /ssc. spitzer. caltech.edu/tools / specpet / ) . 
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which corresponds to roughly 90kpc at z ~ 2, and so 
each target is fully enclosed in the observed slit and is 
thus treated as a point source. With shared coverage of 
the first LL order (LL1: 19.5-38.0/im) and the second LL 
order (LL2: 14.0-21.3/im) we aimed to cover rest-frame 
emission longwards of 6 /xm to search for PAH emission 
and for any silicate absorption at 9.7 /im. With a redshift 
range of z = 0.65-3.2 for the full sample, the LL1 IRS 
module alone provides the relevant coverage for SMGs 
with z > 2.1; while the LL2 module is also required for 
SMGs with z < 2.1. We observed each target at two 
different nod positions for a total integration time of 1-2 
hours with each module. The data were obtained be- 
tween August 2005 and June 2006. 

We reduced the data using the Spitzer IRS S14 
pipeline. 10 To ensure homogeneity throughout the sam- 
ple, we reprocessed the raw data fro m our earlier re- 
sults ( Mc nendez-Delmestre et al"1l2007f) . which had been 
processed by the S13 pipeline. We performed addi- 
tional cleaning of the 2D spectra using the SSC utility 
IRSCLEAN 11 to remove rogue pixels, by creating a mask 
for each Astronomical Observing Request (AOR). We use 
the rogue pixel mask associated with the observing cy- 
cle as a base and then identify additional deviant pixels 
in a superdroop file resulting from the averaging of all 
droop files within a single AOR. We relied on differencing 
between the two nod positions to subtract the residual 
background. To do this, we median-combined all cycles 
for science images of a given SMG with the same module 
order (LL1 or LL2) and nod position. For each mod- 
ule order, we then subtracted the median of all science 
images at one nod position from that of the other nod 
position. 

For each LL order, we used the SPitzer IRS Cus- 
tom Extraction 12 (spice) software to extract flux- 
calibrated ID spectra from each 2D sky-subtracted 
median-combined nod. We used the optimal extraction 
mode included in spice. Optimal extraction entails tak- 
ing a weighted average of profile-normalized flux at each 
wavelength t o increase the S/N in IRS obse rvations of 
faint sources (jHorne 1986; N arron et al.ll2006r ). Flux cali- 
bration is based on observations of a standard star during 
the same observing cycle as the science observations and 
has a 10% uncertainty (see IRS Data Handbook 13 for de- 
tails). The SPlCE-extracted ID spectra for the different 
nods within one LL order were then averaged together. 
We combined the orders using iraf to produce a final 
spectrum for each target, excluding the noisy edges at 
^obs ;$ 14.2/im and X i, s > 35/im. The resulting spectra 
are shown in Fig. [2] 

3. ANALYSIS 

The spectra in Fig. [2] show a range of emission and 
absorption features. One of the key advantages of our 
survey over previous studies is the availability of spec- 
troscopic redshifts for the entire sample from C05. With 
spectroscopic redshifts in-hand, we have a-priori knowl- 
edge of the observed wavelengths where emission and ab- 
sorption features are expected. We therefore identified 

10 http:/ /ssc. spitzer. caltech.edu/irs/dh/ 

11 http:/ /ssc. spitzer. caltech.edu/archanaly/contributed/irsclean 

12 http:/ /ssc. spitzer. caltech.edu/postbcd/spice. html 

13 http:/ /ssc. spitzer. caltech.edu/irs/dh/ 



PAH features, silicate absorption centered at ~ 9.7^m 
and a number of possible forbidden narrow emission lines, 
such as [Aril] (6.99^m) and [Siv] (10.51/W). We see that 
a large fraction of our sample shows PAH emission fea- 
tures and we mark these and the other spectral features 
on Fig. [2] 

We fit our spectra to derive quantitative estimates of 
the strength of the emission and absorption features and 
determine their precise wavelengths. However, given our 
modest integration times (Table [1]), we reach a signal to 
noise {S/N) of only 1-5 in the continuum emission at 
24/im 14 . At this S/N we fo und that intricate spectral 
fitting programs (e.g. pahfit; iSmith et al.ll2007f l are not 
appropriate. We thus adopted a simple and straight- 
forward fitting approach which we describe below. 

Measurements of fluxes and equivalent widths (EW) 
of features in the spectra require that we first determine 
the continuum level for each spectrum before measuring 
the strength of each individual emission and absorption 
feature. The resulting EW and flux measurements are 
extremely sensitive to the fitting method selected, par- 
ticularly in the determination of the continuum level. A 
number of alternative fitting methods have been used 
by different authors, making a fair comparison of PAH 
strengths i n different samples a rather subjective and del- 
icate issue (jSaiina et al.l 120071 ). In particular, the preva- 
lence of interactive fitting methods has further compli- 
cated the reproducibility of mid-IR properties between 
samples. Hence in this work we have adopted a method 
of continuum and PAH fitting that ensures complete re- 
producibility. Moreover, in order to make fair compar- 
isons between the SMGs in our sample and other low- and 
high-redshift sources, we have re-measured all the mid- 
IR features in all of these samples using our method. 
We have applied our method to local ULIRG s from 
the B right IRAS samp l e stud ied by lArmus et all (|2006l 
l2007h and iDesai et at] (|2007h . the low-re dshift nuclear 
starbu rst-dominated galaxies presented bv lBrandl et al.l 
(l2006f) and t he mid-IR s elected high-r edshift sources of 
ISaiina et all (|2007l ) and iFarrah et aD (|2008f ). We note 
that within the context of PAH luminosities and EWs, we 
only include into our discussion those objects for which 
we could distinguish the PAH features. As we mentioned 
earlier, we have also included in our compari sons the 13 
SMGs in the GOODS-North field studied bv lPope et all 
(2008), nine of which have optical redshifts measured by 
C05 and thus naturally complement our sample. 

The thermal mid-IR continuum, arising either from 
starburst or AGN activity can be described by a power- 
law, S v ~ i/~ Qmi r, where the power-law index cumir in- 
dicates the steepness of the spectrum. Flatter mid-IR 
continua with low power-law indices (q?mir ^ 0.5) are 
usually associated with unobscured AGNs, where emis- 
sion from the hot dust close to the AGN is readily ob- 
served. Steeper continua (omir ^ 0.5) are on the other 
hand associated with starburst activity, but can also be 
found in dust-obscured AGNs. Hence we determined the 
power-law continuum index «mir for the spectra in our 
sample to: (1) characterize the continuum slope and in- 
vestigate the underlying energetics, and (2) to define the 
continuum level on which PAH emission features are su- 

14 S/N was determined within a A-window of AA ~ 1.4/xm on 
the individual nod-subtracted 2D spectra 
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Fig. 2. — Mid-IR spectra of all 24 SMGs in our sample redshifted to the rest-frame according to their rest-frame optical redshifts (C05) or 
their PAH-based redshift in the case of SMM J105151.69, SMM J105158.02, SMM J105227.58, SMM J221733.91 and SMM J221806.77 (see 
£|4.1l for details). The spectra are smoothed to the instrument resolution (AA ~ 0.3/im at A ~ 24/xm). The continuum fit for each spectrum 
is denoted by a short-dashed line and the shading represents the flux uncertainty, as given by SPICE. We do not detect SMMJ030231.81 
at z = 1.316. The locations of PAH features at 6.2, 7.7, 8.6, 11.3 and 12.7 (im are indicated by the dashed lines and these features are 
visible in a large fraction of our sample. The position of tentative narrow line emission from [Aril] (6.99/jm) and/or [Siv] (10.51/xm) are 
shown by dotted lines in SMM J105158.02, SMM J105238.19, SMM J123711.98, SMM J123721.87, SMM J163639.01, SMM J221733.91 and 
SMM J221735.15. 

continuum for each spectrum are shown in Fig.O 

Having determined the continuum, we search for the 
presence of PAH emission features at 6.2, 7.7, 8.6, 11.3 
and 12.7/im. Since our objects cover a range in red- 
shifts of z ~ 0.65-3.2, the wavelength coverage varies 
from ~ 9-22/im for the lowest redshift SMG to ~5-9/im 
for that with the highest redshift. Out of all detected 
SMGs in our sample, >90% cover of the 7.7/zm region of 
the spectrum, thus providing a large sample of 21 SMGs 
to investigate the spectral features within this impor- 
tant wavelength region. We define a wavelength window 
around each spectral feature and use the iraf task SPLOT 
to fit the PAH features with individual Gaussians, where 
the center and FWHM are left as free parameters within 
the fitting window. We create a cursor file specifying the 



perimposed. 

To determine «mir we use the x 2 minimization IDL 
routine LINFIT to obtain the best linear fit to each SMG 
spectrum in log(Sv)-log(i/) space. We first obtain the 
best linear fit to each full spectrum and calculate the rms 
in the (data — fit) residual spectrum. Then, we identify 
those wavelength regions in the original spectrum where 
the data points lie more than 1 x the rms away from this 
first linear fit. In this way, we automatically identify 
the wavelength regions where emission and absorption 
features are present. We then perform a second round of 
linear fitting, where we exclude these wavelength regions 
from the fit. We adopt the results of this second fit as 
our final continuum, with the influence of emission and 
absorption features minimized. The resulting fits to the 
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fitting windows and continuum level for each spectrum 
to automate the fitting process and avoid the lack of re- 
producibleness inherent to other interactive fitting tasks. 
The wavelength regions used to fit the 6.2, 7.7, 8.6, 11.3 
and 12.7//m PAH features were approximately 5.90-6.60, 
7.10-8.3, 8.30-8.90, 11.00-11.60, 12.40-13.00 /xm, respec- 
tively. We adopt the SPLOT output uncertainties in the 
resulting PAH EWs and integrated fluxes, which account 
for the noise level present in the spectra. These results 
are presented in Table O 

We take advantage of the PAH features in our spectra 
to determine mid-IR redshifts for the SMGs in our sam- 
ple. For this we derive a redshift associated to each PAH 
feature from its observed central wavelength. The final 
PAH-based redshift for spectra displaying multiple PAH 
features corresponds to the weighted mean of the red- 
shifts derived from the individual PAH features, where 
the weights are provided by the PAH fitting uncertain- 
ties. We present these results in §4.11 

Finally we quantify the strength of the silicate ab- 
sorption. We define the optical depth as rg.7u m = 



Continued. 

\og[S con t/ S b s ], where S b s is the flux level observed at 
9.7/im and S con t, the interpolated flux of continuum in 
the absence of any absorption. We measure the latter 
from our power law continuum fit. In ! )4.4I we present 
the resulting values for the SMGs in our sample that 
provide sufficient wavelength coverage for a reliable de- 
termination of the continuum. 

3.1. Composite Spectra 

In the interest of obtaining a representative spectrum 
for the sample as a whole, with a higher S/N to en- 
able identification of fainter underlying features, we also 
median-combined the individual rest-frame spectra of the 
SMGs in Fig. [2] into a composite spectrum. The full col- 
lection of individual spectra in our sample cover the rest- 
frame wavelength range A ~ 5-22 /im. To avoid being 
overly biased towards the mid-IR properties of SMGs 
within a narrower redshift range, we only include in 
the composite spectrum wavelengths covered in at least 
five SMGs. The individual spectra in Fig. [5] also dis- 
play a significant spread in rest-frame continuum level 
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Fig. 3. — (Left:) The median-combined rest-frame composite spectrum of 22 SMGs detected in our sample, with the shaded area 
representing the 1-cr sample stan dard deviation. We assume optical redshifts (C05), except for the five SMGs in which PAH features 
suggest an alternate redshift (see !|4.H . The composite spectrum is dominated by strong PAH emission at 6.2, 7.7, 8.6, 1 1.3 a nd 12.7 (im, 
with an underlying red continuum. All the spectra are normalized by their flux at A = 8.5 /im in the rest-frame (see t!3 - 1 1 for details). 
(Right:) Composite spectra for subsets of SMGs in thr ee separate redshift bins. These demonstrate that at A > lOfim our composite 
spectrum is dominated by sources at z < 1.6. See ^4.5. II for discussion. 



with Sg^m ~ 0.03-0.2 mjy and so we normalize them by 
their flux at a rest-frame wavelength of 8.5/im. Only the 
spectrum of SMMJ221733.12 does not cover this wave- 
length and is thus excluded from the composite. The fi- 
nal composite spectrum thus includes 22 detected SMGs 
and is shown in Fig. [3] The sample standard devia- 
tion of the composite spectrum increases towards longer 
wavelengths, in part reflecting the fact that fewer in- 
dividual spectra contribute at X re st ^ 10/im. To bet- 
ter demonstrate the redshift ranges which contribute to 
the different wavelength regions in the composite (and 
their variations) we divide our sample into three differ- 
ent redshift bins and construct composite spectra of low- 
, intermediate- and high-rcdshift SMGs (see Fig. [3]) and 
discuss these in §4.5.11 

We also take advantage of the range in properties 
covered by our sample (see Fig. []} to investigate mid- 
IR spectral properties of different SMG sub-populations. 
We follow the method described above to construct com- 
posite spectra of subsamples of SMGs based on 24/im- 
and radio brightness to search for trends within our 
sample. However, composite spectra based on apparent 
brightness mix sources with different intrinsic luminosi- 
ties at different redshifts. We attempt to minimize this 
bias by first defining luminosity-complete subsamples of 
SMGs, then restrict the SMGs included in each compos- 
ite to those that fall within a redshift range where we 
are close to complete in terms of mid-IR and radio lumi- 
nosity. At the median redshift of our sample, <z>^ 2, 
the observed 24/nm flux corresponds to redshifted rest- 



frame 8/mi emission. Therefore, we define luminosity- 
complete subsamples of SMGs by using the rest-frame 
luminosities at 8/im and at 1.4 GHz. 15 . To build the 
24/im-composites we restrict the parent sample to SMGs 
with L 8Mm > 10 32 erg s -1 . For the radio composite, we 
include SMGs with L 1AGHz > 2 x 10 24 ergs" 1 Hz" 1 . As 
a result, low-redshift SMGs with the lowest mid-IR and 
radio luminosities in the full sample, even though they 
have the highest apparent fluxes, are excluded from the 
24fim- and radio-brightness composites. Results for the 
various composite spectra from §3.1 are presented in Ta- 
ble [31 We discuss these composites in H4.5.2I 

4. RESULTS 

Out of the 24 SMGs in our sample shown in Fig. [21 we 
detect PAH emission and continuum in 15 SMGs, PAH 
emission and no continuum in four, continuum but no 
PAH emission in another four and neither continuum nor 
PAH emission in only a single SMG (SMM J030231.81 at 
z = 1.316). Hence, our sample contains 19 SMGs with 
PAH emission and 19 with detectable continuum emis- 
sion. In addition, marginal detections of [Aril] (6.99/im) 
and [Siv] (10.51/im) narrow-line emission were found in 
several individual SMGs (see Fig. 

4.1. PAH-based Redshifts 

16 hi.iGHz = ^D 2 L S 1A GHz{l + z)- (Q + 1) WHz-\ assuming 
a = —0.8, the av erage spectral index for star-forming galaxies (e.g., 
lYun et al.ll200lD 
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Fig. 4. — A comparison of the PAH-derived redshifts with the 
rest-frame UV redshifts for the SMGs in our sample (C05). PAH- 
based redshifts, with a typical uncertainty Az = 0.03 agree well 
with the optical redshifts determined by C05, except for five out of 
the 23 detected SMGs, namely SMM J105151.69, SMM J105227.58, 
SMM J105158.02, SMM J221733.91 and SMM J221806.77. These 
SMGs have multiple PAH features that suggest a redshift signifi- 
cantly different from the optical one. Four out of these five SMGs 
have mid-IR spectra displaying one particularly prominent feature 
that we associate with the 7.7fim PAH line. The presence of ad- 
ditional fainter features confirm the validity of this assumption. 



We have rest-frame UV redshifts for all the SMGs 
in our sample from optical spectroscopy using the 
Low- Resolution Imaging Spectrograph on Keck II (C05). 
These spectroscopic redshifts allowed us to identify con- 
fidently faint emission and absorption features in the 
mid-IR spectra of individual SMGs. Nevertheless, with 
19 SMGs displaying prominent PAH emission, we can 
determine a mid-IR redshift for these SMGs based on 
the location of the redshifted PAH features and compare 
these to the rest-frame UV redshifts. PAH-based red- 
shifts (zpah) for our sample are shown in Tableland are 
compared to the redshifts from C05 in Fig. 0] We note 
that PAH-derived redshifts may have large uncertainties 
(Az ~ 0.01 — 0.05) arising from the intrinsically large 
width of PAH features (FWHM> lO^ms- 1 ). Redshifts 
based on rest-frame UV features are thus more precisely 
determined. However, we note that because the optical 
emission in these systems does not trace the bolomet- 
ric emission well, there are opportunities for these rest- 
frame UV redshifts to be erroneous if the optical spectro- 
scopic IDs are spatially offset from the radio counterpart. 
Equally, the large IRS aperture may result in mid-IR 
bright sources unrelated to the SMG to dominate the 
spectrum. 

The comparison in Fig. |4] demonstrates that the IRS 
spectra of the majority of the SMGs have features con- 
sistent with the rest-frame UV redshifts. However, there 
are five SMGs where the features suggest an alternate 
redshift (see TableEJ that is significantly different from 
the rest-frame UV redshift. For the five SMGs with dis- 
crepant redshifts, the reasons for the disagreement are in 
all cases due to ambiguities and line misidentifications re- 
sulting from low S/N features in the UV spectra, and so 
we discuss these individually. We note that these correc- 



tions to the redshifts have not significantly changed the 
interpretation of these sources. 16 We also note that one 
of the new PAH redshifts, SMM J105158.02, strengthens 
the z ~ 2.7 SMG association discussed in iBlain et al.l 
(|200l . 

In the case of SMM J221733.91, an initial indication 
of a redshift around z ~ 2.5 from absorption features in 
a low S/N spectrum was apparently verified by a weak 
Ha line at z = 2.555 (jSwinbank et al.ll200l ). However, 
a deeper UV-spectrum, doubling the S/N, revealed an 
unambiguous z = 2.865 from multiple absorption fea- 
tures, consistent with the z = 2.89 found from the IRS 
spectrum. 

Similarly with SMM J105227.58, the original noisy UV 
absorption line spectrum was somewhat ambiguous but 
suggested z — 2.142. Subsequent deeper spectroscopy 
has revealed this source to be z — 2.470 from weak Lya 
in emission as well as Civl549 absorption and emission, 
consistent with the PAH redshift within fitting uncer- 
tainties. 

SMMJ221806.77 was identified in C05 at z = 3.623 
from a single weak emission line detection in a relatively 
sky-free spectral region and corroborated by a red con- 
tinuum dropping out in the _B-band. The PAH-based 
redshift z = 1.91 suggests that the continuum of this 
source is heavily reddened, and that the rest-frame UV 
emission feature is either noise or a superposed line de- 
tection from a different galaxy. 

SMMJ105158.02 was identified at z ~ 2.2 first in the 
UV through various low S/N absorption features. An 
apparent Ha detection cons istent with z = 2.23 9 was 
proposed as the true redshift (jSwinbank et al.ll2004h . Re- 
analysis of more recent, deeper UV-spectra does reveal 
a suite of absorption features consistent with z — 2.694 
(Sin 1260, Oi/Sinl303, Cnl335). 

Finally, for SMM J105151.69, an apparent emission line 
proposed as [On] 3727 is not reproduced in later Keck 
spectra taken on this galaxy at a different position an- 
gle, and it has now been identified as an emission line 
galaxy slightly offset from the radio position. The PAH- 
based redshift of z = 1.62 is consistent with UV-spectra 
[Cm] 1909 absorption features and Mgn 2800 emission 
features. 

4.2. PAH Luminosities and PAH ratios 

The 7.7/xm PAH feature is the most prominent in the 
mid-IR spectra of star-forming galaxies and is thus par- 
ticularly useful in studying star formation in faint sources 
at high redshift. We derive PAH luminosities from the 
integrated line fluxes presented in Tableland plot these 
in Fig. [5] We find the bulk of 7.7/im PAH luminosities 
for SMGs in the range ~ 10 44 -10 45 ergs _1 and hence 
that, as expected from their far-IR luminosities, SMGs 
have PA H luminosities that e xceed those of low-redshift 
ULIRG s (lArmus et al.l I2007T) 17 and local nuclear star- 
bursts (jBrandl et al.l l2006f ) by factors of up to 10 and 
1000, respectively. These PAH luminosities indicate that 
~ 0.5-5% of the bolometric luminosity in SMGs escapes 

16 Changes in redshift at the 10-20% level at z ~ 2 do not result 
in significant changes to the luminosities or temperatures beyond 
the associated calibration errors. 

17 We include NGC 6240, Mrkl014, UGC5101, Arp220, 
FSC 05189-2524, FS C 12112-0305 , Mrk 273, FSC 14348-1447, 
FSC 22491-1808 from lArmus et al.l [12001 ). 
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Fig. 5. — The 7.7/im-PAH luminosities as a function of redshift 
for the SMGs in our sample. In addition we plot the SMGs from 
Pope ct al. (2008). For comparison we also plot apert ure- corrected 
low-r edshift nuclear starburst-dominated galaxies (Brandl^et^l 
2006|), loca l ULIRGs from the bright IRAS sample (lArmus et~ 
2006I. I2007P and high-red shift mid-IR-selected sources (Saii na et al 
20071 ; IFarrah et al.ll200g|) . We see that SMGs are among the most 
PAH-luminous objects in the universe, with PAH features up to 
1000 X and 10 X more luminous than local starbursts and ULIRGs, 
respectively. Considering that PAH emission is associated to strong 
star-forming activity, this suggests that the SMG population com- 
prises some of t he most vigorous star formation events ever ob- 
served. See i|4.2l for discussion. 

in the form of PAH emission. 

From Fig. [5] we can see that the range in PAH lumi- 
nosities for SMGs covers that occupi ed b y the mid-IR 
selecte d sources in lSaiina et al.l f2007f) and IFarrah et al.l 
(|2008f) . The sample from ISaiina et ail (|2007l ). based 
on mid-IR colors and bright 24/im fluxes (S^vm > 0.9 
mjy) , has been shown to comprise mostly c omposite ob- 
jects, dominated by AGN activity (~ 75%; ISaiina et al.1 
120071) . Roughly 25% of their sample display strong PAH 
features, which explains the large PAH luminosities dis- 
played in Fig. ranging from values comparable to 
low-redshift ULIRGs to those shared with SMGs. Even 
though these objects have similar PAH luminosities to 
that of SMGs, their PAH EWs are on average lower, 
with SMGs displaying a median 6.2/im PAH EW three 
times larger and a median 7.7/xm PAH EW larger by 
a factor of ~ 1.5. Probing slightly lo wer 24/im fluxes 
(S24fim ^ 0.5 mJvL lFarrah et al.1(|2008t) also select bright 
mid-IR sources at high redshift, but with an additional 
criterium of stellar-dominated near-IR SEDs. These 
objects concentrate around z ~ 1.7 and Fig. [5] shows 
that they display similar PAH luminosities to SMGs 
and to the most PAH- luminous objects in the sample 
fromlSaiina et ai] (12007!) . The mid-IR select ed objects in 
IFarrah et al.l (|2008( ), in contrast to those in lSaiina et al.l 
(|2007l ). have PAH EWs similar to SMGs and more than 
a factor of two larger than t hose of the mid-IR bright ob- 
jects in ISajina et al.1 (|2007l ). indicating a reduced AGN 
contribution to t heir mid-IR emissi on compared to the 
latter sample (see lFarrah et a l. 2008). We discuss further 
the comparison of SMGs and other galaxy populations 
in EJ5J 

Different PAH emission features arise from dis- 



tinct bending /vibra tional modes of the PAH molecules 
(jDraine fc Li 2007). The different modes can be en- 
hanced relative to each other as a result of varying PAH 
ionization state and PAH size distribution. Exposure to 
a more energetic radiation field may strip PAH molecules 
of peripheral H-atoms. A decrease in the fraction of 
neutral PAHs would diminish the strength of the 8.6, 
11.3 and 12.7/mi PAH features, which are produced by 
the in-plane (8.6/iin) and out-of-plane (11.3, 12.7/zm) 
C-H bending mode. As a result, the 6.2 and 7.7/im 
PAH features, which are thought to be produced by C- 
C skeleton vibration (jDraine fell 121)071 ). would appear 
to dominate. Within this context, the relative strength 
of PAH features can be used to p robe the energetics o f 
the underlying radiation field fe.g.. lGalriano et al.ll2008h . 
The size distribution of PAH molecules may also affect 
the observed relative strength of features: smaller PAHs 
Nc 5, 10 3 ) tend to emit m ore strongly at 6.2 and 7.7/im 
Allamandola et al"1ll989f ). In the hottest regions close 
to an AGN, sublimation of these smaller PAH molecules 
may result in a supression of short- wavelength PAH emis- 
sion. Dust extinction has also been shown to play an 
important role in explaining the variation in the rela- 
tive strength of PAH features. In particular, strong ex- 
tinct ion can increase the 7.7 /6.2 and 7.7/8.6 PAH ra- 
tios (jRigopoulou et al.lll999f ) due to an ice absorption 
featu re at ~ 6/im (whi ch reduces the 6.2/im PAH emis- 
sion, ISpoon et al.ll2002l ) and broad silicate absorption at 
~ 9.7/im, where the former reduces the strength of 6.2/im 
PAH feature and the latter that at 8.6/im. 

We consider the flux ratios in the 6.2, 7.7 and 11.3/im 
PAH features in our sample to investigate the energetics 
and extinction in the star-forming environments within 
SMGs. We plot these in Fig.[6]and find that SMGs have 
median 7.7/11.3 and 7.7/6.2 PAH ratios of 3.8 ± 1.6 and 
2.5 ± 1.4, respectively. The median 7.7/11.3 PAH ra- 
tio is similar to that of low-redshift nuclear starbursts, 
S7.7/S11.3 = 3.7. However, the distribution in values 
for the 7.7/11.3 PAH ratio in SMGs is very broad, ex- 
tending over a similar range to that of the high- r edshift 
24/i m-selected sources studied by ISaiina et al.l (|2007l ) 
and IFarrah et al.1 ((2008) . Such divers ity in PAH ratios 
is also present in local U LIRGs (e.g., iRigopoulou et al.1 
Il999t IPeeters et al"1l2004l ). indicating that diverse radia- 
tion fields are found amongst both low- and high-redshift 
IR-luminous sources. We note that the need for measur- 
able 11.3/im emission means that this comparison only 
includes SMGs at z < 1.6 (Fig. EJ. 

In Fig. [6] we illustrate that the bulk of SMGs have 
7.7/6.2 PAH ratios lower than the median value found for 
local nuclear starbursts (St.t/Sq.2 = 2.5 for SMGs vs. 5.3 
for the starbursts). This suggests that an intrinsic differ- 
ence in either radiation field or extinction may exist be- 
tween SMGs and low-redshift nuclear starbursts, which 
result in a suppressed 6.2/im PA H feature and a relativel y 
prominent 7.7/im PAH feature (jRigopoulou et al.ll!999f ) . 
As we discuss later, we believe that reduced extinction 
may be responsible for this behaviour in SMGs. 

4.3. Continuum Slopes 

The continuum power-law index oimir, characterizes 
the steepness of the mid-IR continuum and requires ad- 
equate wavelength coverage to be measured reliably. In 
Table [4] we present the «mir values for the SMGs whose 
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Fig. 6.— Distribution of 7.7/6.2 and 7.7/11.3 PAH ratios within 
the SMG population. Objects for which we have only measur- 
able 7.7/6.2 PAH ratios are plotted at SV.7/S11.3 < and those 
for which we only have 7.7/11.3 measurements are plotted at 
S7.7/S6.2 < 0. We also include PA H ratios measured for a sample 
of low-redshift nuclear star bursts jjBrandl ct al. 2006 ), and high- 
rcdshift mid-IR selected sources (Sajina ct al. 2007; Farrah ct al. 
1200a) . SMGs share similar 7.7/11.3 PAH ratios with these pop- 
ulations, but their distribution in 7.7/6.2 PAH ratios is skewed 
towards values that are in general lower than that of low-redshift 
starbursts. We suggest that this difference in 7.7/6.2 PAH ratio 
may arise due to a difference in extinction (see i|4.2l and £|5.2I for 
discussion). 

spectra cover a wavelength range that extends longwards 
of the silicate absorption and that we thus consider re- 
liable. We find that the SMGs in our sample have a 
median mid-IR power-law index of < otyaR >— 1-05, 
typical of the continuum found for star-forming regions. 
However, the individual spectra span a range of power- 
law indices cumir ~ —0.8-2.5. Considering the uncer- 
tainties in the individual fits, with a median average 
of < Ackmir >^ 0.3, this range in power-law indices 
suggests a real diversity in the mid-IR continuum slopes 
within the sample. Power-law indices in ~ 60% of the 
SMGs in our sample are consistent with those found in 
star forming regions («mir ^ 0.5), while the remaining 
~ 40% have cumir < 0. 5, more charact eristic of AGN- 
dominated regions (e.g. JDeo et"aHl2007l) . 

For a subset of our sample we can also characterize 
the continuum shape from the Se^m/ S 12 color, which 
is measured at wavelengths devoid of substantial contam- 
ination from PAH emission. We indicate these results in 
Table H] At 6/zm, the continuum traces thermal emis- 
sion from hot (T^t ~ 500 K) dust, while at 12/im, the 
continuum is dominated by dust heated to lower tem- 
peratures (T dust ~ 250 K). We find that ~ 60% of the 
SMGs for which we can determine mid-IR colors show 
red continua, with Se^m/ S 'n^m 0.3 and Omir ^ 0.5. 
These resemble the colors of low-redshift star-forming 
gal axies and th e nucl ear starburst galaxies presented 
bv iBrandl efafl (|2006f ). with S 6fMm / S 12fim ~ 0.18-0.53. 
The remaining 40% of this subsample have bluer con- 
tinua {S 6 ^ m /S 12 ^ m ) ~ 0.8-1.1 and a M m < 0.5), more 
typical of the high -redshift mid-IR selected sources of 
ISaiina et all (|2007l ) with S 6lim /Si^ m > 0.3. 

Taken together these results suggests that the mid- 



IR continuum emission in ~ 60% of the SMGs in our 
sample has the spectral characteristics associated with 
star forming systems locally. The remaining third of the 
SMGs show more AGN-like mid-IR continua, although 
these frequently also display strong PAH emission. 

4.4. Obscuration as Measured by Silicate Absorption 

One of the clearest indicators of absorption in the 
mid-IR arises from intervening amorphous silicate dust 
grains. Silicate absorption from the stretching of the 
Si-0 bonds is centered around 9.7/xm and is the main 
absorption feature within the wavelength range consid- 
ered in this work. It is measurable in SMGs at z < 2.0, 
corresponding to roughly half of our sample. 

We use the strength of the silicate-absorption feature 
as a measure of the obscuration along the line of sight to 
the warm continuum emission, Tg.7^ m . We note that this 
definition of optical depth assumes that the obscuring 
material lies in a foreground screen, which is likely to be 
too naive a model for the complex mix of absorption and 
emission in SMGs. 

We calculate Tg^^m for the individual SMGs, but stress 
that since we only detect faint continuum with S/N< 
3 for a number of SMGs in our sample, the values we 
derive for the strength of the silicate absorption are likely 
lower limits (see Table [2]). We find a range of optical 
depths limits Tg.7 Mm ~ 0.05-1.2, with a median value 
of Tg.7^ TO = 0.31. 18 For the composite spectrum, which 
benefits from an improved S/N, we find an optical depth 
of r 9 . 7Ajm = 0.44 ± 0.06. 

In Fig. [7] we present optical depths (rg.7^ m ) as a func- 
tion of the 7.7/im PAH luminosity relative to the hot- 
dust continuum, £7.7^/ Xs^. We compare SMGs to 
low-redshift starbursts and ULIRGs, and to high-redshift 
mid-IR selected sources. We stress that the IRS obser- 
vations of the local ULIRGs contain 80-90% of the star- 
formation activity in these systems and hence the differ- 
ences in apparent properties are not due to comparing 
galaxy-integrated measurements for the SMGs to mea- 
surements of , just th e nuclear component in local systems 
(jArmus et al.ll2007f ). 

The PAH-to-continuum luminosity ratio closely resem- 
bles a measure of EW as it probes the strength of the 
PAH feature with respect to the continuum at neighbor- 
ing wavelengths. The distribution in PAH-to-continuum 
values for SMGs is on average si milar to that of mid - 
IR selected high-redshift sources of iFarrah et al.1 (|2008h . 
but higher than th at of the mid-IR bright sources of 
ISaiina et all (|2007t ). This is in agreement with differ- 
ences in their PAH EWs as discussed in §4.21 suggesting 
an increased AGN contribution to the mid-IR in the lat- 
ter sample. 

Fig. [7] also shows that at low redshifts, starburst galax- 
ies with the highest PAH-to-continuum ratios display 
the highest optical depths, suggesting that the most ac- 
tive starbursts are also the most obscured, with optical 
depths c lose to that of loc al ULIRGs, < Tg. 7pm >= 0.7 
(see also lDesai et alJl2007t ). SMGs lie off this trend, dis- 
playing a distribution in PAH-to-continuum values sim- 
ilar to low-redshift starbursts (Lr.r/Le ~ 0.1-0.8), but 
maintaining on average rg. 7Mm limits that are typically 



18 I ncluding the SMGs with redshifts from C05 from Po pe et al.l 
||2008| ) extends this range to rg.7 Mm ~ 1.4. 
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Fig. 7. — A plot of the 7.7pm PAH luminosity relative to 
the 6-pvn continuum luminosity as a function of silicate optical 
depth. We include optical depth measurements only for those 
SMGs that provide sufficient wavelength coverage for reasonable 
conti nuum estimates: 10 SMGs from our sample and seven SMGs 
from P ope et al.l 1120081) . Optical dep ths for local nuclear starbursts 
HBrandl et al.112009) , local ULIR Gs IIArmus et alll20(}7l) and high- 
redshift 24^m-selected sources (Saiina et al. 2007; Farra rTet al.l 
I2008D are also shown for comparison. Accurate measurements of 
optical depths require careful determination of the unobscured con- 
tinuum. We thus quote the values for high-redshift sources display- 
ing low S/N continuum detection at A = 9.7pm as lower limits. Op- 
tical depths measurements for objects for which we do not have con- 
tinuum measurements at Sa^m are plotted at Log(Lr^ / Lg) > 1; 
upper limits for Log(L-j .7 / 'Lq) in SMGs with no rg.7 Mm measure- 
ments are shown at rg.7 Mm < —0.1. The error bar displays the 
median uncertainty in the SMG PAH-to-continuum ratios. The 
histogram at the top shows the distributions for SMGs and the 
two low-redshift samples. SMGs display rg.7 Mm limits lower than 
local ULIRGs and the most obscured starbursts, but have propor- 
tionally stronger 7.7fim PAH emission. 

less than those seen in low-redshift ULIRGs and in the 
most obscured local nuclear starbursts. This difference 
in distribution suggests that SMGs have lower extinction 
along the line of sight to their warm continuum emission 
sources with respect to these low-redshift samples. Such 
a difference in distribution of silicate optical depths with 
respect to the loc al samples is also se en in the mid-IR 
selected objects of iFarrah et al.l (|2008l ). 

4.5. Composite Spectra 
4.5.1. The SMG Spectrum 

With clear PAH features in > 80% of the SMGs in our 
sample, the composite spectrum in Fig. [3] displays promi- 
nent PAH emission. Taking advantage of the wide red- 
shift range covered by our sample, we explore the com- 
posite spectra of SMGs in three different redshift bins: a 
low-redshift bin {z < 1.6), an intermediate-redshift bin 
(z = 1.6-2.4) and a high-redshift bin (z > 2.4); see Fig.H 
The SMGs in the low-redshift bin dominate the compos- 
ite spectrum at A > 10/im. Therefore, we can track the 
evolution of the mid-IR spectral properties of SMGs at 
A < 10/im using our sample. 

All three redshift composites in Fig. [3] display strong 
7.7/im PAH emission and the 8.6/zm feature also re- 
mains prominent. With increasing redshift we observe 



a broadening of the 7.7/im PAH feature and a weaken- 
ing of the 6.2/mi PAH feature. The increase in scatter 
around the median composite spectrum indicates either 
that a more diverse range in spectral properties is present 
at the higher redshifts covered by our sample or poorer 
S/N. However, we note that the three composites dis- 
play an increase in continuum level with redshift. Within 
the context discussed in H4.2[ the weakening of the ob- 
served 6.2/im-PAH emission feature may be an indication 
that SMGs at higher redshifts suffer from increased ex- 
tinction, hampering detection of PAH emission features 
at the shorter wavelengths, or from a hardening of the 
radiation field the PAH molecules encounter (see also 
IWatabe et ai]|2009fl . 

4.5.2. Investigating Potential Sample Selection Biases 

To investigate how representative our conclusions are 
for the whole SMG population we need to investigate 
the variation of mid-IR properties of SMGs with both 
apparent 24/im and radio flux. Considering that one of 
our selection criteria for building the SMG sample re- 
lied on choosing objects with higher (estimated) 24/im 
fluxes, we probe for differences between 24/im-bright 
and 24/im-faint SMGs (in the observed frame) to un- 
derstand how representative our sample is of the SMG 
parent population with the faintest 24-/mi counterparts. 
Using the SMGs within the Lg^m-complete sample de- 
scribed in H3.ll we construct mid-IR composite spectra 
based on 24/mi-brightness: SMGs with SW^m ^ 0.4 mJy 
comprise the 24/im-bright composite, while those with 
SiA^m ^5 0.4 mJy make up the 24/im-faint composite. 
We plot these composites in Fig. [5] 

We see strong 7.7/im PAH emission in the 24/im-faint 
composite with EW7.r^ m = 2.48 ± 0.35, while a greater 
contribution from hot dust continuum in the 24/im- 
bright composite results in a significantly lower 7.7/^m 
PAH EW, EW 7 .7 M m = 0.99 ± 0.14. Thus it appears that 
the strong PAH emission we see in our SMG sample at 
z ~ 2 is not likely to be a result of our 24/im-criterium 
in defining our initial target sample: the composite spec- 
trum of SMGs with faint 24/im-fluxes also displays strong 
PAH features, even stronger than those in the 24/im- 
bright sample. The latter point even suggests that the 
prevalence of sources with strong continuum emission 
(and hence the contribution from AGN or highly ob- 
scured starbursts) is likely to be higher in our sample 
than in the general SMG population. 

We use the range in radio fluxes covered by our sam- 
ple (see Fig. [T]) to investigate possible intrinsic trends 
in mid-IR properties associated with radio-brightness. 
We seek to address potential differences between the 
C05 sample of radio-identified SMGs and the < 30% 
of 

•SssOjum ^ 5 mJy SMGs with radio-counterparts below 
the sensitivity of current radio surveys. It is possible that 
SMGs with no detectable radio counterparts lie at sig- 
nificantly higher redsh ifts than the radio-id entified pop- 
ulation (z > 3, C05; lYounger et al.ll2007ft . We cannot 
test this suggestion with our sample as we have targeted 
SMGs with known redshifts. However, it has also been 
suggested that the radio-undetected SMGs lie at simi- 
lar redshifts to the dominant radio-detected population, 
but that they have somewhat lower IR luminosities and 
slightly colder characteristic dust temperatures, a com- 
bination which would result in similar observed 850-/^m 
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Fig. 8. — (Left:) The composite spectra of the seven 24/im-bright SMGs within the Lg Mm -complete sample (bottom) and of the seven 
24itm-faint SMGs (top). (Right:) Composite spectra of the five radio-bright SMGs within the Li.4G^f z -complete sample (bottom) and of 
the five radio-faint SMGs (top). We normalize the spectra of the individual SMGs that comprise the radio composites at A = 7.8/^m. The 
24/im-faint composite is significantly more PAH-dominated than the 24/im-bright composite, while the latter displays a stronger continuum. 
In the case of the radio-composites, the radio-bright sample appears to display a broader 7.7/xm PAH feature over a slightly more prominent 
continuum. However, characterization of the radio-composite spectra does not reveal any significant difference within our measurement 
uncertainties (see Tabled- 



fluxes (|Chapman et alJ l2004). Following this argument, 
our radio-detected sample may therefore be biased to- 
wards systems with hotter characteristic dust tempera- 
tures, either due to more compact star forming regions 
or the presence of an AGN. To investigate this possi- 
bility, we divide our sample into radio-bright and radio- 
faint SMG subsets. We first define a luminosity-complete 
subsample of SMGs based on their rest-frame radio (1.4 
GHz) luminosity and use their apparent radio fluxes to 
identify the radio-bright and radio-faint SMGs for which 
we can then construct composite spectra. These are 
shown in Fig. [H] 

The radio-faint SMG composite has a narrower 7.7^m 
PAH feature that emerges sharply from the continuum, 
while the radio-bright composite displays a broader, less 
distinct PAH feature. Considering the line-to-continuum 
parameter traditionally used to gauge the strength of 
PAH features (Z/c; e.g., iGenzel et al.lll998[ ). the radio- 
faint composite displays a higher 7.7/zm-(Z/c) relative 
to the radio-bright composite. A lower 7 .7 )j,m-(l / c) for 
the radio-bright composite may suggest a larger obscured 
AGN contribution. However, the EWs of the 7.7/zm PAH 
feature in these composites are not significantly differ- 
ent within the measurement uncertainties (see Table|3]). 
Therefore, overall our observations do not reveal any sig- 
nificant difference in mid-IR properties based on radio- 
luminosity and we conclude that our radio-detected SMG 
sample is likely to be representative of the whole SMG 
population at z < 3.5. 



4.5.3. Starburst and AGN -dominated Composites 

The SMGs in our sample have been classified prior to 
this work as either AGN- or star-formation-dominated 
or starburst) systems using deep X-ray observations 



Alexander et alj 120051 ) , rest-frame optical spectroscopy 
Swinbank et al.ll2004 iTakata et al"1l2006ft and/or rest- 



frame UV spectroscopy (C05). Table 0] provides a sum- 
mary of these classifications. We construct composites 
for an AGN and a star-forming class based on these clas- 
sifications to investigate differences in mid-IR spectral 
properties and plot these in Fig. [9l Since we do not take 
mid-IR properties into account to build these composites, 
the differences are not the result of any mid-IR selection 
bias. 

We find that the composite spectra of both the AGN- 
and the starburst class display prominent PAH features. 
However, the PAH emission in the AGN composite lies 
on top of a stronger continuum, with a lower 7.7^m PAH 
EW compared with that in the star-forming compos- 
ite (EW 7 . 7t , m ^AGN = 1-25 ± 0.18^m vs. EW 7 . 7Mm ,s B 
2.0!) ± 0.30/im). The AGN composite spectrum also has 
an emission feature at A = 10.51/zm, which is much less 
obvious in the star-forming-composite. This feature may 
correspond to [Siv] narrow-line emission, which is com- 
monly associate d with hot starburst activity or the pres- 
ence of an AGN (|Spoon et al.ll200l lOrie et al.l[2006l) . We 
note that all of the individual SMGs which went into this 
AGN-composite and whose spectra cover the A ~ 10.5/im 
wavelength region show what could correspond to [Siv] 
emission (see Fig. [2]) . Since the SMGs that make up this 
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Fig. 9. — (Left:) The c omposite spectra of eight SMGs in the sample with AGN signatures either in the X-ray (Alexander ct al. 2005), 
optical (C05) or near-IR (Swinbank ct al. 2004) (bottom) and of the 10 SMGs with mid-IR spectra that suggest the presence of an AGN, 
either EW7.7 Mm < 1 or Qmir < 0.5 (top). (Right:) Composite spectra of 16 SMGs in our sample with no AGN signature either in the 
X-ray, optical or near-IR (bottom) and of the 13 SMGs with mid-IR spectra dominated by strong star formation activity, with EW7.7^ m > 1 
and steep mid-IR continuum, ctMiR > 0.5. We see that the AGN-composites display an enhanced continuum relative to the star-forming 
counterparts. We also note that the AGN-composite with AGN signatures at wavelengths ot her tha n the mid-IR (bottom left panel) 
displays a prominent emission feature, likely corresponding to [SIV] emission from an AGN. See ^14. 5.31 for details. 



AGN composite have AGN signatures in the rest-frame 
UV and optical, it is likely that our line of sight has direct 
access into the broad line region. Thus it is also possible 
that this feature corresponds instead to emission from 
solid silicates very near a central AGN. 

Fig. [5] also shows the composite spectra for AGN- 
and starburst- subsets based on their mid-IR classifi- 
cations. We classify SMGs with EWy^™ > 1 and 
steep mid-IR continua (Se^m/ Si2fj,m < 0.3) as star- 
formation-dominated systems. SMGs with EW7.7^ m < 1 
or flat mid-IR continua (Senm/iZfim > 0.3 or q;mir < 
0.5) are classified as 'intermediate'. SMGs displaying 
a continuum-dominated spectrum, either with no sign 
of PAH features or with a continuum-absorption fea- 
ture such as SMM J163650.43 (see Fig.©, we classify as 
'AGN' to indicate a more significant AGN contribution 
(see Table HJ). As expected from following these mid-IR 
classification criteria, the AGN composite in Fig. [5] dis- 
plays a 7.7/im feature with a lower EW than the star- 
burst composite (EWr.r^m^AGN — 1.57 ± 0.22/mi vs. 
EW7,7p m ,sfl — 1-87 ± 0.26^im). We observe a similar 
trend for the other PAH EWs, in particular for the 6.2/im 
PAH feature, which appears in Fig.[9]to have been heavily 
diluted by the enhanced continuum in the AGN compos- 
ite. 

The AGN- and starburst classes, either constructed 
from mid-IR criteria or from classifications based on ob- 
servations at other wavelengths, appear to differ mainly 
in the mid-IR continuum level. The AGN composites 
display an enhanced continuum with respect to the star- 



burst composites, resulting in lower 7.7/im PAH EWs. 
We note that the AGN composites also display promi- 
nent and luminous PAH features. Even though the spec- 
tra of individual SMGs in our sample display 7.7/im PAH 
features spanning an order of magnitude in luminos- 
ity range, the AGN composites have PAH luminosities 
within a factor of ~ 1.5 of those in the starburst com- 
posites. We return to discuss these results in the context 
of the PAH dilution and PAH destruction scenario in |j5l 

4.6. The AGN Contribution to the Bolometric 
Luminosity in SMGs 

In iMenendez-Delmestre et all (|2007l ) we presented the 
analysis of five SMGs that we also include in this pa- 
per. We found that the composite spectrum is well fit 
by a starburst-like template (represented by M82) with 
an additional continuum component, likely arising from 
an AGN. We found that the mid-IR spectral features of 
SMGs are more comparable to the starburst M82, than 
to other well-studied local ULIRGs, such as Arp220 and 
Mrk 1014. We now consider our complete sample of 24 
SMGs and compare the composite SMG spectrum to the 
co mposite of nearby nu clear starburst-dominated galax- 
ies (jBrandl et al.ll2006l : see Fig. [10]). By using this star- 
burst template rather than M82, we seek to represent 
to a better degree the diversity in continuum slopes and 
PAH strengths in starburst-dominated galaxies in the lo- 
cal Universe. 

We find that the mid-IR composite spectrum of the 
SMGs in our sample is reasonably well fit at A < 10/im 
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Fig. 10. — The composite spectrum of 22 SMGs detected in 
our sample as shown in Fig. [3] We fit this composite with the 
mid-IR comp osite spectrum of 13 starburst-dominatcd low-redshift 
galaxies from Brandl ct al. (2006), scaled down to match the SMG 
composite spectrum at 7.7 fim. A power-law component with mid- 
IR spectral slope, omir ~ 1.97, best represents the continuum in 
the composite SMG that is unaccounted for by the starburst tem- 
plate. This continuum may arise either from an obscured AGN 
or optically-thick emission around star-forming regions. See Sec- 
tion l4.6l for details. 

by the local starburst template, but we find that an ad- 
ditional power-law continuum is necessary to fully de- 
scribe the SMG composite spectrum: at the longer wave- 
lengths, A > 10/mi, the starburst template fails to match 
the continuum level (see Fig.[T0|. To quantify this addi- 
tional continuum that is unaccounted for by the starburst 
template we fit a power-law to the (SMG composite — 
starburst-template) residuals. We find that the resid- 
uals are best fit by a power-law F v ~ 1 /-° ! mir ) where 
ttMiR = 1-97 ± 0.22. This additional warm continuum 
may be associated to obscured AGN activity or to an 
additional starburst continuum component. We derive 
an upper limit to the contribution from an AGN to the 
SMG mid-IR spectrum by making the conservative as- 
sumption that this additional continuum arises uniquely 
from a dust-enshrouded AGN. 

Hard X-rays emission provides one of the most direct 
routes to estimate the luminosities of AGN. However, in 
the presence of high colum n densities, such as esti mated 
in SMGs (N H > 10 23 cm" 2 : [Alexander et aHl20f)5h . hard 
X-ray photons may be completely absorbed. In these 
cases an alternate probe for AGN emission is the mid-IR 
emission from hot dust near the AGN. The mid-IR thus 
provi des a complementary insight into a d eeply obscured 
AGN (ISturm et al.ll2006 t (Lutz et al.l l2004h. In particular, 
iKrabbe et al.l ( 2001[ ) show that the 10/im-flux is tightly 



correlated with the X-ray emission, S^-iokeV' both in 
Seyfert and in starburst galaxies, with markedly distinct 
mid-IR slopes. We use the derived correlation for Seyfert 
galaxies to estimate the X-ray luminosity of the AGN 
component in the SMG composite spectrum and hence 
gauge the AGN contribution to the bolometric luminos- 
ity of SMGs. 

From the 10/im continuum flux and the Krabbe et al. 
relation we estimate an X-ray luminosity, Lx ~ 8 x 
10 43 erg s _1 , for the SMGs in our sample at <z>^ 2.0. 



Encouragingly, this estimate is in close agreement with 
the absorption-corrected rest-frame 0.5-8 keV X-ray lu- 
minosities found for the SMGs in the Chandra Deep Field 
(| Alexander et al.l I2005D . < L x >= 5 x 10 43 erg s _1 . 
To assess the contribution of the AGN to the energet- 
ics of the SMGs we first estimate the far-IR luminosi- 
ties of the SMGs using the observed radio-far-IR cor- 
relation for SMGs (jKovacs et al.ll2006l ). yielding an av- 
erage far-IR luminosity of £f ar tr ~ 5 x 10 45 e r g s" 1 . 
We then use the approach of Alexander et al.l (|2005l ) 
and lMenendez-Delmestre et al.l ( 2007t ). and adopt the X- 



ray-to-FIR luminosit y ratio typically found for quasars 
(Lx /^far-iR ~ 0.05; I Alexander et al.ll2005h to estimate 
that the AGN in SMGs contribute ~ 32% of the SMG's 
far-IR luminosities or Lf ar _iri ~ 1.6 x 10 45 erg s _1 . 

We emphasize that this AGN contribution estimate is a 
strong upper limit for the SMG population. Considering 
that the enhanced continuum in the mid-IR bright com- 
posite suggests that our selection criteria may bias us to- 
wards an increased AGN contribution within our sample 
(see M4.5.2|) . the AGN contribution estimate based on the 
SMG composite is likely biased high. It is also very likely 
that some fraction of the red continuum unaccounted for 
by the local starburst template arises from dust emis- 
sion heated by optically-thick starbursts in SMGs, par- 
ticularly considering the high mid-IR power-law index, 
o:mir = 1.97 ± 0.22, in the range associated with regions 
dominated by star-formation. Therefore, we quote 32% 
as a firm upper limit on the AGN contribution to the 
bolometric luminosity in typical SMGs. 

With an SMG composite dominated by PAH features, 
it is clear that the bolometric luminosity of SMGs is in 
general dominated by star-formation. However, the in- 
dividual mid-IR SMG spectra display modest variation 
in AGN contribution and in some cases can reveal the 
presence of an AGN that remains invisible at other wave- 
lengths (see Table|4j. In particular, there is evidence that 
SMMJ123553.26 is a Compton-thick AGN. This SMG 
is undetected in the X-ray but is bright and dominated 
by continuum emission in the mid-IR, thus displaying 
no mid-IR signatures of star formation (see Fig.[2|). Ex- 
trapolating from its mid-IR spectrum, we estimate a hot 
dust luminosity 19 , LQ^ m ~ 2.0 x 10 45 erg s" 1 . Fo llowing 
the approach described bv I Alexander et all (|2008l ) we use 
the 0.5-2 keV observed luminosity of SMM J123553.26 
(|Alexander et al.l [20 05) to derive a rest-frame 1.6-24.8 
keV luminosity and convert it to an intrinsic X-ray lumi- 
nosity, i2-iokov < 2 x 10 42 erg s -1 , assuming an X-ray 
spectral index, T = 1.4. On the basis of the X-ray-6/im 
luminosity of this SMG, we find that it lies significantly 
below the X-ray-6/zm luminosity relationship found for 
nearby AGN, suggesting that t his source hosts a lum i- 
nous Compton-thick AGN fe.g.. I Alexander et aLl r2008). 
However, we note that this SMG lies very close to a sec- 
ond radio source (< 5"), which complicates the deblcnd- 
ing of th e individual flux con tributions in MIPS 24/im- 
lmagmg (jHainline et al.ll2009f l. We detect a single con- 
tinuum trace in the IRS 2D spectrum that, even though 
it coincides with the position of SMM J123553.26, likely 
corresponds to the blending of the mid-IR emission of 
both sources. Therefore, it remains a possibility that the 
measured mid-IR flux of this source does not correspond 
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to the X-ray source. 

We compare the IR luminosity associated with the 
AGN component in the SMG composite to that of ob- 
scured extremelY_iumincjas quasars at z = 1.3 — 3 stud- 
ied by Pollctt a'et al.l ((2008). Since the latter have been 
shown to have < 20% of their bolometric luminosity aris- 
ing from starburst activity, it is reasonable to assume 
that the bulk of their IR luminosity arises from AGN 
activity. We find that the IR emission contributed by 
the AGN in SMGs is roughly 10% the IR luminosity 
of these luminous quasars. Since the bolometric lumi- 
nosities in these two samples are similar, the respective 
AGN contributions to the total IR luminosity conform 
with the evidence that these two populations are dom- 
inated by different power sources. Furthermore, this 
demonstrates that systems dominated by star-formation, 
such as SMGs, may attain as high luminosities as AGN- 
dominated systems at high redshift. 

4.7. Comparison to SMGs in Other Samples 

With a s ample of 23 detected SMGs in thi s work, plus 
nine from (Valiantc et al. 20071) and 13 from (|Pope et al.l 
2008), we have statistics that allow us to uncover the 
range in mid-IR properties present within the SMG pop- 
ulation. The majority of SMGs display prominent PAH 
features atop an underlying weak power-law continuum, 
potentially arising from an AGN. These spectra suggest 
starburst-dominated activity, with varying AGN contri- 
butions. 

We find only one clear example of a continuum- 
dominated source with a prominent feature centered at 
A ~ 8/xm from unabsorbed contin uum, SMM J163650.43 
(|Menendez-Delmestre et alj|2007j ). Three other SMGs 
in our sample display rather featureless continuum- 
dominated spectra (SMM J105238.30, SMM J123553.26, 
and SMMJ221737. 39). In the sample o f ni ne detected 
SMGs presented bv lValiante et al.l (|2007l ) and lLutz et all 
(|2005h . SMMJ02399-0136 displays PAH features over- 
laying a particularly strong continuum and is thus clas- 
sified as a source powered by equal contri butions of star 
form ation and a Compton-thick AGN (jValiante et alJ 
[2001 . Out of the 13 S MGs in the GOOD S-North field, 
SMMJ123600.15 (CI; iPope et alJ 120081) also displays 
PAH features on a steeply rising continuum and is clas- 
sified as an SMG with a 44% AGN contribution to the 
bolometric luminosity. 

We note that 25% of the SMGs in our sample lie at 
rcdshifts z < 1.2, which is a redshift range not covered 
by the other SMG mid-IR studies. This extension to 
lower-redshifts provides us with an insight into the SMG 
emission at the longer mid-IR wavelengths, A > 11.5/um. 
In Fig.Hniwe can see that it is at these longer wavelengths 
where the local starburst template becomes increasingly 
insufficient to account for the SMG mid-IR emission and 
where the additional continuum component becomes in- 
creasingly significant. Coverage of this wavelength region 
within our sample is thus crucial to compare the com- 
posite SMG spectrum with potential local analogs and 
hence to assess the AGN contribution, which we mea- 
sure to be <32% of the bolometric luminosity. Coverage 
out to these longer wavelengths is also critical for better 
constraining the mid-IR continuum slope, giving addi- 
tional insight to the nature of the underlying power en- 
gine based on the measured steepness of the mid-IR spec- 
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Fig. 11. — PAH luminosities as a funct ion of IR luminosi ty for 
the SMGs in our sample and for those in Pope et al. (2008|). We 
also include apertu re-corrected PAH luminosities of low-redshift 
nuclear starbursts (Brandl et al. 2006). The line shows the best 
linear fit to both low-redshift and high-redshift objects (see Eq.fl}. 
At low redshifts the 7.7/xm luminosity correlates well with the IR 
luminosity. At higher redshifts, SMGs appear to mantain such 
correlation. 

trum. Improved continuum constraints allows us in turn 
to quantify the strength of the 9.7-/im silicate absorption 
and investigate the distribution of silicate dust in SMGs, 
compared to other galaxy populations (see ^5.2[) . 

5. DISCUSSION 

We compare the mid-IR properties of the SMGs in our 
sample with those of starb urst-dominated low-redshift 
galaxies ([Brandl et al.l |2006[). select local ULIRGs from 
the Bright IRAS sample (jArmus et al.l 120061 l2007f) and 
the mid-IR s e lected sources of ISaiina et alJ (|2007l) and 
iFarrah et all ([20081 ) . We use PAH luminosities, PAH 
relative intensities, the steepness of the continuum and 
the strength of silicate absorption to compare the level 
of star-forming activity, the energetics in radiation envi- 
ronments and the distribution of cold and warm dust in 
these galaxies. 

SMGs have among the most luminous PAH emission 
in the universe, displaying levels of activity which are 
unseen in local samples (Fig.[5]). This result reflects the 
significant increase in star-formation activity in the high- 
redshift SMG population which underpins the thousand- 
fold increase in the space density of ULIRGs out to z ~ 2 
(C05). The similarly large PAH luminositie s displayed by 
the m id-IR selected high-redshift objects in lFarrah et al.l 
(2008) also put in evidence the marked increase in PAH 
luminosities with redshift compared to local ULIRGs. 

The presence of large PAH luminosities (£7.7^™ > 
10 44 erg sec -1 in SMGs with lower PAH EWs and/or 
flatter mid-IR continua, such as SMM J105238.30 and 
SMMJ123549.44, and in SMGs with AGN signatures 
at other wavelengths, such as SMM J123711.98, suggests 
that even in the presence of an AGN, PAH molecules are 
abundant in these galaxies. This is also suggested by the 
large PAH luminosities displaye d by the AGN - domi nated 
composite objects discussed in ISaiina et al.l (|2007t ) and 
shown in Fig. O It has been suggested, under the PAH 
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destruction scenario, that highly energetic photons from 
an AGN overcome the PAH binding energy and reduce 
the populat ion of PAH molecules and the observed PAH 
luminosity (|Voitj|1992T ). However, rather than a decrease 
in PAH luminosity, we find that SMGs with flatter mid- 
IR continua and AGN signatures at other wavelengths 
tend to have lower PAH EWs. Similar results have been 
rec ently discussed in the context of low-redshift ULIRGs 
bv lDesai et al.l (|2007ft . where a decrease in EWs for the 
6.2 and 11.3/im PAH features is observed with increasing 
24/im rest-frame luminosity. These results suggest that 
an increased AGN contribution results in an enhanced 
dust continuum that dilutes the prominence of PAH fea- 
tures with respect to the local continuum, rather than in 
a destruction of PAH molecules. Our study allows us to 
extend the validity of the PAH- dilution scenario out to 
the higher redshifts of SMGs. 

5.1. PAH Emission as a Star Formation Rate Indicator 

In the past, authors have explored the use of the mid - 
IR light as a measure of galaxy SFRs. lElbaz et alJ (|2002l ) 
claim that the flux measured in the 6.75/xm ISO band is 
related to the total IR luminosity, which in the absence 
of an AG N is taken to be a good tracer of s tar- fo rming 
activity (iKennicuttl 1 1998ft. Roussel et al.l (|2001l ) and 
iForster Schreiber et al.l ( 20041 ) showed that there exists 



a correlation between ISO 6.75//m light and Ha luminos- 
ity in local spirals and starburst galaxi es. With the ad - 
vent of Spitzer, 8/mi IRAC resul ts from lWu et al.l (|2005l ) 
and I Alonso-Herrero et a l. (2006) further confirmed that 
the 8/im band flux is closely associated with the star- 
formation taking place in low-redshift (z < 0.2) luminous 
infrared galaxies with Lir < 10 12 L . 

More recently, several authors have exploited the 
6.2/im PAH emission feature as a SFR indicator, based 
on the observed correlati on between the 6.2/xm PAH and 
total IR luminosity (e.g., iBrandl et al.ll2006l : iPone et all 
2008). The wavelength coverage within our study gives 
us an insight to the 7.7/im PAH emission feature of 22 
SMGs. In Fig.[TT] we show the 7.7/rni PAH luminosi- 
ties as a function of IR luminosities for SMGs and low- 
redshift starburst galaxies. We see that L 7 . 7tim and L/jj 
for low-redshift starbursts are well correlated and that 
the SMGs extend this local correlation out to higher red- 
shifts, where we have assumed that the bulk of the total 
IR luminosity can be associated to star formation (see 
q4.6[) . We find that the best fit for the complete sample 
of objects at low- and high-redshifts is: 



io 5 (L 8 _iooo M m) = (0.99±0.05)xZo 5 ( J L 7 . 7Aim ) + (2.45±2.3). 

w 

The spectral insight to the mid-IR emission of galax- 
ies allows us to disentangle the PAH emission from the 
AGN contribution to the underlying continuum, which 
remains otherwise intermixed in mid-IR broad band 
photometry. In the presence of an AGN, the contin- 
uum contribution may be substantial and even domi- 
nate the emission at these mid-IR bands (see for exam- 
ple the continuum-dominated SMGs SMM J123553.26, 
SMMJ163650.43 and SMM J221737.39 in Fig.[2|). Tak- 
ing the IR luminosity as a proxy of star-formation activ- 
ity (|Kenni cutt 1998]), we use the relation between PAH 
luminosity and IR luminosity to characterize the 7.7/im 



PAH luminosity as a SFR indicator. Based on this, we 
find that the 7.7/xm PAH luminosity can be used as a 
SFR indicator following the following form: 

log(SFR) = -40.9 + 0.99 x log{L 7 . 7 ^ m ), (2) 

where L-jj^ m and SFR are in units of erg s _1 and 
M Q yr _1 , respectively. The bulk of SMGs have PAH- 
derived SFRs ~ few xlO 2 - 1 3 M^ yr" 1 , consistent 
with the Xray-derived SFRs by Alexander et al.l (|2005l ) 
for SMGs in the CDF-North. We obtain a median value 
of < SFR > 7 7 ».m~ 800 Mq yr " 1 for SMGs, including 
the objects in lPope et all (|2008ft . 

5.2. Dust Distribution in SMGs 

Our observed silicate optical depths correspond to a 
median visual extin ction of Ay ~ 6, assuming Ay — 
(18.5 ± 2) x T 9 . 7fim (|Draindl2003ft . This is significantly 
higher than the values estimated in previous near-IR 
studies, Ay ~ 1 — 4, based on the Balmer decrement 
dTakata et al.ll2006ft and broad-band optical /near-IR col- 
ors dSmail et al.ll2004tlSwinbank et alJ200l . This differ- 



ence demonstrates that these optically-based techniques 
are only capable of probing the less-obscured regions 
of the SMGs and hence return lower values for Ay. 
The mid-IR gives us a supplementary insight to the 
star-formation taking place within the highly optically- 
obscured regions of SMGs and allows us to quantify the 
additional extinction at rest-frame optical wavelengths 
unaccounted for by the Balmer decrement. The differ- 
ence in values for Ay derived from the mid-IR and the 
optical/near-IR indicates that Ha-derived SFRs are (on 
average) suppressed by a factor of > 100 when uncor- 
rected for internal extinction. 

The range in silicate optical depths we find in 
SMGs, T9.jfj.rn ~ 0.05-1.2, falls below the typical val- 
ues reported fo r both lo cal ULI RGs and AG Ns (< 
T 9 .7 M m >~ 1. 5: IShi et all 12006] lHao et alJ 12007ft. Re- 
cent work by iDesai et al.l ( 20071 ) measured the strength 
of silicate-absorption features in a large sample of low- 
redshift ULIRGs from the Bright IRAS Sample. They 
find that cold ULIRGs, with steeper mid-IR spectra 
Szsftm/ Smpim S 0-2, have deeper silicate-absorption fea- 
tures than warm ULIRGs, with S^s/mi/S'eo^m 0.2. We 
find that SMGs have 79 7 /im -values close to that of warm 
ULIRGs, such as 72L45FSC05189-2524, and to the me- 
dian value found for star-forming low-redshift galaxies 
(< T- g .7 Mm >= 0.47 ± 0.38; see Fig. 0). 

The strength of the silicate feature can provide in- 
sight into the distribution of dust along the line of sight. 
The strong silicate absorption measured for cold low- 
redshift ULIRGs can be explained by a small, deeply 
embedded source with large amounts of obs curing dust 
along the line of sight (|Levenson et al.ll2007ft . Shallower 
silicate-absorption features and the presence of silicate- 
emission, such as found in t he high-redsh i ft m id-IR- 
selected sources presented by ISaiina et alJ (|2007ft and 
warm low-redshift ULIRGs suggest that the mid-IR has 
a direct view to the hot thermal continuum source, pos- 
sibly through a clumpy obscuring medium. 

The observed 7.7/6.2 PAH ratios of SMGs (Fig.EJ) are 
on average lower than those seen in local nuclear star- 
bursts. This ratio is potentially sensitive to ionization 
and reddening, with softer radiation fields and lower ob- 
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scuration leading to lower 7.7/6.2 ratios. When taken 
with the differences in silicate absorption between SMGs 
and local ULIRGs, both these results can be explained 
by a single cause: lower extinction towards the mid-IR 
line and continuum emitting regions within SMGs, com- 
pared to either local nuclear starbursts or ULIRGs. The 
activity in these latter populations are concentrated into 
highly compact regions and the lower extinction in the 
SMGs suggests the mid-IR visible star formation in these 
galaxies is likely to occur in a more extended component, 
more similar to the activity in the disks of "normal" star- 
forming galaxies in the local Universe. Similar evidence 
for a more physically extended distribution in starburst- 
dominated g alaxies at z ~ 1.7 has been reported by 
iFarrah et all (j2008l ). 

Evidence for extended star formation within SMGs has 
also been found at other wavelengths. High-resolution 
MERLIN/ VL A observations of two samples of SMGs 
have revealed radi o morphologies extending out to ~ \" 
or~ 8kpc in size (IChapman et al]l2004t iBiggs fc Ivisonl 
I2008f) . These two studies find that ~ 85% of the com- 
bined sample of 24 SMGs have extended radio emis- 
sion, while only ~ 15% appear to be dominated by an 
unresolved component. Adopting radio emission as a 
proxy for far-IR emission, these results suggest that mas- 
sive star-formation is occurring on ^ 8 kpc scales within 
SMGs. This is in contrast to results for local ULIRGs, 
where the far-IR emission is confin ed to a compact nu- 
clear region of ~ 1-2 kpc in size (jCharmandaris et alJ 
2002). Further evidence for physically extended star 
formation within SMGs has co me from studies of their 
rest-frame optical emission lines. ISwinbank et~aTl (|2006l ) 
mapped Ha emission from eight SMGs with Integral 
Field spectroscopy and found Ha structures extending 
to 1-2" or ~ 8-16 kpc in size. These results differ, how- 
ever, from those of parallel mapping of the molecular 
gas emission in SMGs. Subarcsecond resolution millime- 
ter imaging of CO emission from six SMGs with the 
IRAM Plateau de Bure Interferometer suggests compact 
gas distributions in these galaxies, with sizes < 2 kpc 
(jTacconi et al.ll2006l |2008| ) . These observations typically 
target the high-J CO transitions, which trace warm and 
dense gas, so it remains possible that there are more 
extended reservoirs of cold gas on ^ 2 kpc scales within 
these systems which could fuel their extended star forma- 
tion. Our mid-IR results give further evidence towards a 
more extended distribution of the star-forming regions in 
SMGs, in contrast to the compact nuclear starbursts seen 
in similarly luminous local galaxies. These very extended 
starbursts may reflect the structural properties of the 
progenitors of SMGs, which are more gas-rich and may 
have smaller bulges than those which form local ULIRGs. 

6. SUMMARY AND CONCLUSIONS 

SMGs are an enigmatic population which have been 
proposed to be undergoing tremendous amount of star- 
forming activity, in order to explain their ultra-luminous 
infrared emission. The high stellar masses and SRFs in- 
ferred for these galaxies have prompted the idea that 
SMGs are the likely progenitors of today's most massive 
galaxies. However, the contribution from AGN activity 
to the bolometric luminosity has remained a caveat to 
these results. We present the largest sample of SMGs 
observed in the mid-IR with Spitzer IRS. Our main re- 



sults are: 

1. At the individual level, we find that 80% of SMGs 
in our sample display luminous PAH features, with 
£7.7/™ ~ 10 43 -10 45 ergs -1 . This confirms for a large 
sample of SMGs that this is a population undergoing in- 
tense star-forming activity and that it includes some of 
the most intense star-formation events ever witnessed in 
the Universe. 

2. We find only four cases of SMGs with continuum- 
dominated spectra (17% of our sample), which indi- 
cates that even though SMGs appear to be starburst- 
dominated as an ensemble, some diversity is present 
within the sample. 

3. We quantify the strength of the silicate absorp- 
tion feature at A ~ 9.7/im and, keeping in mind that 
such measurement requires a careful determination of 
the continuum, we report lower limits of Tg.7^ m only for 
those SMGs that provide sufficient wavelength coverage 
for reasonable continuum slope estimates. We find that 
the distribution in Tg.j^m limits for SMGs falls below that 
of local ULIRGs and the most obscured low-redshift nu- 
clear starburst-dominated galaxies. This suggests that 
SMGs have lower dust obscuration to their mid-IR con- 
tinuum emitting regions than local ULIRGs. 

4. Comparison of Sr.r^m/ Su^^m PAH flux ratios sug- 
gests that SMGs host similar radiation environments to 
local star-forming galaxies. However, the SV^rn/ 'Sq^hw. 
PAH ratio is lower in SMGs than in local nuclear star- 
bursts. This can most easily be explained if the extinc- 
tion to these regions is lower in SMGs than in the star- 
burst population. 

5. From the composite SMG spectrum we measure the 
potential AGN contribution as revealed by an enhanced 
hot-dust continuum at \ re st ^ 10/im. We make con- 
servative assumptions regarding the AGN contribution 
and find that the relative contribution of an AGN with 
respect to star-formation activity remains low. With a 
maximum AGN contribution of < 32% to the total lu- 
minous output, SMGs are clearly a galaxy population 
dominated by intense starburst activity. 

6. We conclude that the detailed mid-IR spectral prop- 
erties of SMGs show several differences to local ULIRGs 
or nuclear starbursts. These differences can be most eas- 
ily explained by a difference in the extinction to the mid- 
IR continuum and line emitting regions of these galaxies, 
with the SMGs showing systematically lower extinction. 
Given the compact geometry of the mid-IR emission from 
local ULIRGs, the difference in extinction we see ar- 
gues that the mid-IR emission in SMGs, both continuum 
and PAH-emission, arises in a more extended compo- 
nent (> 1-2 kpc). This supports radio and Ha studies 
which indicate that star formation in SMGs is extended 
on scales of > 1 — 2 kpc, far larger than seen in local 
ULIRGs. 
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Individual SMGs 
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a Optical rcdshifts from C05 PAH-bascd Rcdshifts with a median uncertainty of Az ~ 0.03; sec CTTl for details. C PAH fluxes arc expressed in units of 10~ 15 crgs~ 1 cm 2 . Quoted errors correspond to SPLOT PAH fitting uncertainties 
and flux uncertainties added in quadrature. 6 PAH rcst-f'ramc EWs are expressed in units of fiin Optical depth derived from silicate-absorption feature. Sec 94.4| for details. g Fluxes reported are average fluxes within X obs = 23.7 ± 2.35/im 
in the individual IRS spectra, to reproduce the 24^m-MIPS band wavelength coverage.* 1 When PAH features remain undetected above the noise we quote upper limits corresponding to the integrated flux within the PAH fitting windows 
(see J3} in the error spectrum from SPICE. 1 We detect faint continuum in these SMGs^ PAH properties of low-S/N features arc quoted as upper limits. ^ This SMG displays a mid-IR continuum-absorbed bump centered at ~ 8/nm. We quote 
the integrated flux within this bump as an upper limit to the 7.7fim PAH emission.* Only half of the PAH feature falls within the wavelength coverage. 
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TABLE 3 

Composite Spectra 
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Median rcdshift of SMGs comprising composite spectrum. PAH luminosities arc oxprosscd in units of 10 44 ergs" 1 . Quoted errors correspond to SPLOT 

d p 
PAH fitting uncertainties and flux uncertainties added in quadrature. PAH rest-frame EWs are expressed in units of [fim] . S 24Atm > 0.40 mjy from the 

8fim luminosity-complete subsamplc of SMGs with Lg Mm > 10 32 ergs _lf S 24 ^ m < 0.40mJy from the 8fim luminosity-complete subsamplc of SMGs with 

L 8(lm > 10 32 ergs- lg S 14GHz > 95^Jy from the luminosity-complete subsamplc of SMGs with Lj . 4GHz > 2 X 10 24 ergs" 1 Hz" lh S ± . 4GHz < 95/iJy from 

the luminosity-complete subsamplc of SMGs with L x 4 q Hz > 2 X 10 24 ergs" 1 Hz" 11 Composite of SMGs classified as starburst-dominated in the rest-frame 

optical iSwinbank ct al. 2004; Takata et al. 200G), UV (C05) and/or X-ray 1 Alexander et al. 2005) J Composite of SMGs classified as starburst-dominated in the 

mid-IR, with EW 7 _ 7flTn > 1 and c* MIR > 0.5. See text for details. k Composite of SMGs with AGN signatures in the X-ray [Alexander et al.|[2'OT5T ) and/or in 

the ncar-IR jSwinbank et al J| 200^; I Takata et al.1120061) 1 Composite of SMGs with AGN signatures in the mid-IR: E\V 7 7 ^ Lm < 1 and/or c MIR < 0.5. See text 

for details. 
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TABLE 4 

AGN AND STARBURST (SB) SIGNATURES IN RADIO-IDENTIFIED SMGS 
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a We report mid-IR colors for SMGs with wavelength coverage \ res t ~ 6— from [Alexander et al ] J2005h c from C05 d Classification 



based on mid-IR spectral properties (see fl4.5.3| l. c from Swin bank et al ] (I2004fl f from lTakata et al.l l|2006ft 



